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ABSTRACT 



O 

cn 

^ ' Multicolor [BVRI) light curves of the eclipsing classical nova V Per are 

O ■ presented, and a total of twelve new eclipse timings are measured for the 

system. When combined with previous eclipse timings from the literature, 
these timings yield a revised ephemeris for the times of mid-eclipse given by 
JD0 = 2,447,442.8260(1) + 0.107123474(3) E. The eclipse profiles are analyzed 
■ with a parameter-fitting model that assumes four sources of luminosity: a white 

dwarf primary star, a main-sequence secondary star, a flared accretion disk with 



a rim, and a bright spot at the intersection of the mass-transfer stream and the 
disk periphery. Model parameters include the temperatures of the white dwarf 
(Ti) and the secondary star (T2), the radius (i?d) and temperature (Tj), of the 
disk periphery, the inner disk radius (-Rm) , the disk power- law temperature expo- 
nent (a) and thickness (/ir), and a bright spot temperature enhancement factor 
(Xs)- A matrix of model solutions are computed, covering an extensive range of 
plausible parameter values. The solution matrix is then explored to determine 
the optimum values for the fitting parameters and their associated errors. For 
models that treat the accretion disk as a flat structure without a rim, optimum 
fits require that the disk have a flat temperature profile. Although models with a 
truncated inner disk (i?in >> i?i) result in a steeper temperature profile, steady- 
state models with a temperature profile characterized by T(r) oc r~^^^ are found 
only for models with a significant disk rim. A comparison of the observed bright- 
ness and color at mid-eclipse with the photometric properties of the best-fitting 
model suggests that V Per lies at a distance of ~ 1 kpc. 
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Subject headings: binaries: eclipsing - novae, cataclysmic variables - stars: dwarf 
novae - stars: individual (V Persei) 

1. Introduction 

V Per is an 18th magnitude classical nova system (Nova Persei 1887), which reached 
an estimated magnitude of m ~ 4 — 5 at the peak of its eruption (see McLaughlin 1946). 
More than a century later, the system was discovered to be eclipsing by Shafter & Abbott 
(1989) who established an orbital period of 2.57 hr. V Per is of particular interest because 
its orbital period places the system near the middle of the 2-3 hr "gap" in the orbital 
period distribution of cataclysmic variable stars (Kolb 1996; Ritter 1996). In addition to the 
possibility that V Per could have been formed in the period gap, Shafter & Abbott discussed 
possible evolutionary scenarios that could have brought V Per to its current period under 
the assumption that the system evolved from an initially longer orbital period. 

Eclipsing cataclysmic variables are valuable because they provide the opportunity to 
study the the radiative properties of accretion disks through analysis of the eclipse profiles. 
Generally, properties of the system are determined either through parameter-fitting models 
or through maximum-entropy eclipse mapping techniques (e.g. Horne 1985). In the former 
case, properties of the accretion disk, such as its radial temperature profile, must be specified 
in analytic form, while the latter method allows for an arbitrary disk emissivity distribution, 
and selects the "simplest" (maximum entropy) disk emissivity map. 

Shortly after V Per was discovered to be echpsing. Wood et al. (1992) analyzed the 
original monochromatic photometric data from Shafter & Abbott (1989) using the maximum 
entropy eclipse-mapping technique. The principal result of their study was that the accretion 
disk in V Per had a relatively "fiat" radial temperature profile that was considerably less 
steep than the T(r) oc r~^^* profile predicted for a steady-state disk (e.g. Shakura & Sunyaev 
1973). Despite the inconsistency with the expected steady-state relation, the fiat temperature 
profile seen in V Per is not unusual. Several eclipse mapping studies of novalike variables and 
dwarf novae in eruption - systems expected to harbor steady- state disks - have also exhibited 
surprisingly flat temperature profiles (e.g. Rutten et al. 1992; Robinson et al. 1999; Biro 
2000). This discrepancy has been explained by Smak (1994) who used an echpse-mapping 
technique to analyze model eclipse profiles from disks with and without rims, demonstrating 
that the spurious temperature profiles found in previous studies of high inclination systems 
{'^ 80°) were a result of a failure to model the disk rim. More recently, Knigge et al. (2000) 
used synthetic eclipse profiles of a rimmed disk from their study of DW UMa to show that 
self-occultation of the disk by its rim, in addition to the contribution of the rim to the total 
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light, causes the eclipse profiles to be more V-shaped than in the case of a fiat disk. They 
argued following Smak (1994), that if these light curves were then modeled assuming the 
disk is flat (as in eclipse mapping), an erroneously flat temperature gradient would result. 
In the case of V Per, Wood et al. (1992) found that a model with a truncated inner disk, 
such as what one might expect if the white dwarf primary were strongly magnetic, had a 
temperature profile that was considerably steeper, similar to that expected for a steady-state 
disk. 

Tabic 1 summarizes the properties of the 15 cataclysmic variables with measured periods 
between 2.25 hr and 2.75 hr. Two of the stars arc classified as classical novae (one of which is 
V Per), eight as magnetic systems (i.e. DQ Her or AM Her stars), two as novalike variables, 
two as SU UMa dwarf novae, and one as both a nova and a DQ Her star. Given that eight 
out of the 15 systems near the middle of the period gap, are in fact magnetic (i.e. either 
DQ Her or AM Her systems), it would not be surprising if Wood et al. (1992) were correct, 
and V Per turned out to be a magnetic system as well. In order to pursue this question 
further, and to explore how the addition of a disk rim would affect the model results, we 
have obtained multicolor eclipse light curves of V Per and analyzed the data with a new 
parameter-fitting eclipse model. In this paper, we present the results of that study. 



2. Observations 

Observations of V Per were carried out in 2002 October and December using the 2048^ 
CCD imager on the Steward Observatory Bok 2.3 m refiector. A summary of observations 
is presented in Table 2. To increase the time sampling efficiency, only a small subsection 
of the full array was read out. The subsection was centered on V Per and included several 
nearby stars that were used as comparison objects for differential photometry. Typical time 
resolutions were in the range of 40-60 sec. 

The data were processed in a standard fashion (bias subtraction and flat-flelding) using 
IRAF.^ The individual images were subsequently aligned to a common coordinate system 
and aperture photometry was performed through a 6" aperture on V Per and two nearby 
comparison stars, one much brighter than V Per and the other of approximately equal 
brightness. Differential light curves of V Per were constructed with respect to the brighter 
comparison star and the stability of the atmospheric conditions was monitored using the 



^IRAF (Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by AURA, Inc., under cooperative agreement with the National Science 
Foundation. 
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star of comparable brightness. In all cases, the monitor star showed no long-term drift, and 
had a point-to-point scatter of less than 0.02 differential magnitudes. 

To calibrate our observations, several standard stars were observed on each night. Zero 
points were determined on each night using mean extinction values appropriate to the Kitt 
Peak site. On several nights, standards were observed over a sufficient range of airmass for 
a full solution including the extinction coefficient in each band. In all cases, we found ex- 
tinction coefficients consistent with the mean values. All subsequently reported magnitudes 
were computed using nightly determined zero points and the mean extinction. The two com- 
parison stars in the V Per field were calibrated by calibrating every individual frame of our 
observations on each night and computing the mean. The comparison star used to calibrate 
the V Per hght curves is located at RA=2:01:54.3, DEC=+56:44:38.7 (2000.0) (~ 35" north 
of V Per), and is characterized by S = 16.90 ± 0.03, V = 16.11 ± 0.03, R = 15.60 ± 0.03, 
and / = 15.09 ± 0.03. The twelve cahbrated eclipse light curves - three each in B, V, R, 
and / - are shown in Figure 1. V Per is characterized by an out-of-eclipse magnitude of 
V ~ 18.1 ±0.1 and B — V c:^ 0.6± 0.2. The luminosity of the system appears to be relatively 
stable over the ~ 2 month time span covered by the observations, with the exception of the 
data obtained on 29 Oct 2002, when V Per appeared to be ~ 0.4 mag fainter than usual. 
Mean magnitudes and colors of V Per, both in and out of eclipse, are given in Table 3. 

3. Updated Ephemeris 

Eclipse timings for V Per extend back more than 15 years to the original eclipses dis- 
covered in the work of Shatter & Abbott (1988). Table 4 summarizes all available timings of 
mid-eclipse from Shatter & Abbott (1989), Dahm (1997), Katysheva et al. (2002), and from 
the present study. Following earlier studies, times of mid-eclipse for the twelve eclipses in 
the present study have been determined by fitting a parabola to the lower half of the eclipse 
data. A hnear least-squares fit of the mid-echpse times from Table 4 yields the following 
ephemeris for V Per: 

Tmid-eciipse = JD© 2, 447, 442.8260(1) + 0.107123474(3) E. (1) 

Residuals of the individual eclipse timings with respect to this ephemeris are also given 
in Table 4, and are plotted as a function of cycle number in Figure 2. The period of V Per 
appears to be remarkably stable, with no evidence for any significant period change over the 
~ 15 years that the eclipses of V Per have been monitored. 
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4. Eclipse Model 



The V Per data have been analyzed using an improved version of the echpse modehng 
program described in Shafter et al. (2000, Paper I), and in Shatter & Holland (2003, Pa- 
per II). The model considers four principal sources of light from the system: the white dwarf, 
the secondary star, the accretion disk, and the bright spot where the inter-star mass transfer 
stream impacts the periphery of the disk. The model light curve flux at a given orbital phase 
is computed by summing the contribution from the secondary star and the unocculted re- 
gions of the remaining three sources of radiation. The fluxes from all four sources of radiation 
have been corrected via the use of a linear limb-darkening law with coefficients, fix, given in 
van Hamme (1993). The program operates on the normalized light curve intensities, so in 
effect, we fit only the shapes and depths of the eclipse profiles. The program does, however, 
compute a model B — V color from the assumed blackbody fluxes using the transformations 
given in Matthews & Sandage (1963), which can then be compared with the observed B — V 
color. 

The major improvement in our model is in the treatment of the accretion disk, which is 
no longer confined to the orbital plane. The disk thickness perpendicular to the orbital plane 
now varies linearly with radius, resulting in a rim at the disk periphery. The flaring angle 
of the disk and the rim (half) height are parameterized by a fltting parameter, hj- — h/Rd, 
where h is the disk height at its edge perpendicular to the orbital plane, and is the disk 
radius. As before, the disk is assumed to radiate like an optically-thick blackbody, with a 
temperature distribution that follows a radial power-law profile given by: 



where is the temperature of the rim at the outer edge of the disk. For an optically thick, 
steady-state disk, a characteristic value of a = 0.75 is expected (Shakura & Sunyaev 1973). 

Another improvement to earlier versions of our model involves the treatment of radiation 
from the secondary star. In the current model, the tidally-distorted shape of the Roche- 
lobe-filling secondary star is fully taken into account, not only when computing the occulted 
regions of the white dwarf, disk, and bright spot, as in earlier versions of the code, but also 
when computing the contribution of the secondary star to the model fiux. The fiux from the 
white dwarf component is computed assuming it to be a spherical blackbody of temperature 
Ti. The radii of the white dwarf and secondary stars are computed from the orbital period 
and mass ratio as described in Paper I. 

One of the significant advantages of the inclusion of a disk rim is that the bright spot 
can be more easily modeled. The presence of the disk rim provides both an isotropic and 
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an anisotropic component to the bright spot radiation, so there is no need to parameterize 
these components separately as in Paper II. The isotropic component of the bright spot 
is produced as in Paper I, by radiation from a region of the disk's surface defined by the 
intersection of a circular area of radius 0.2 (centered on the point of intersection of 
the mass-transfer stream and the disk perimeter) and the accretion disk. The anisotropic 
component is produced by radiation from the corresponding azimuthal region of the disk rim. 
As before, the bright spot is assumed to radiate like a blackbody, the bright spot temperature 
parameterized by a multiplicative factor, Xs, applied to the local disk temperature in the 
bright spot region. 

For a given mass ratio and orbital inclination, the geometry of the eclipse is determined 
by dividing the primary's Roche lobe into a matrix of cells having Cartesian coordinates 
{x, y, z), and then computing an occultation kernel 0{x, y, x, 0) for eclipse by the lobe- filling 
secondary star. The x and y axes lie in the orbital plane with the positive x axis pointing 
toward the secondary star, and the positive y axis oriented 90° counterclockwise as seen from 
the direction of the positive z axis. For a given orbital phase, (p, the occultation kernel takes 
on values of zero or one, depending on whether or not the cell is eclipsed. The monochromatic 
fluxes from the accretion disk at orbital phase are then given by the following summations. 
For the disk surface and disk rim, respectively, we have: 

. x=Rd y=\f^^ 

/a(</') = ^E E ll{x,y)^x^yO{x,y,z^K^f^^\^), (3) 

. . <t>'=-K/2 z=hdRd 

/a(0) = ^ E E Il{x,y) {RaA<f>') cos <f>' AzO{x,y,z,<f>), (4) 

0'=— 7r/2 z=-hrRd 

where d is the distance to the system, i is the orbital inclination, and (f)' = (f) + tan~^(y/x) is 
the phase-modified disk azimuth angle. The Cartesian coordinates of the disk rim are given 
hy X — Ra cos (0' — (f)) and y — sin (0' — 0). 

The disk surface and rim are assumed to radiate as blackbodies, with limb-darkened 
specific intensities given by: 

Il{x,y) = B^{T^S^,^[Rj^x^ + yY) (1 - /^a + /^acos 7), (5) 

y) = BxiT^S^^y) (1 - //A + //A sin i cos 0'), (6) 

where cos 7 = cos (3 cos i — sm. i sin /3 cos 0' is the cosine of the angle between the line-of- 
sight and the normal to the disk surface, (3 — tan~^(/ir) is the half-angle of the disk fiare, 
and //a are the linear limb-darkening coefficients. The parameter 5'x,y modifies the outer 
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disk temperature in the vicinity of the intersection of the mass transfer stream and the disk 
periphery {xs,ys) by the factor Xs ^ follows: 



Xs when i/^a^^a^^j^TT^^-y^ < ^spot 
1 otherwise. 



The fluxes for the secondary star and the white dwarf primary are computed as follows. 
For the secondary star we have: 

^aI*^) = ^ E Bx{T2)[l -l^x + /^Acos C(0)] cos AA, (7) 

where AA is an element of area on the surface of the tidally-distorted secondary star with a 
normal vector inclined at an angle ^ to our line of sight. The white dwarf primary is assumed 
to be a spherical blackbody of temperature Ti and radius Ri. Its integrated flux is given by: 

fli<P) = ^5a(Ti) (1 - /iA/3) g(i?in) TcRl 0(0, 0, 0, 0), (8) 

where the parameter Q applies a correction depending on whether the lower hemisphere of 
the white dwarf is shadowed by the inner edge of the disk, Ria, as follows: 

i±f^ when i?in = Ri 
1 when i?in > Ri- 



Q{Rin) = 

Finally, the normalized light curves, independent of the distance to the system, are given 



by: 



where a normalization phase of 0n = 0.2 has been adopted to represent the out-of-eclipse 
hght level. 



4.1. Input Parameters 

The input parameters required by the model arc summarized in Table 5. They can 
be divided into two general categories: fixed parameters and fitting parameters. The fixed 
parameters are those that are either known, a priori, such as the orbital period, or param- 
eters whose values can be estimated using assumptions implicit in the model. The fitting 
parameters, on the other hand, cannot be specified initially, and are varied during the fitting 
procedure. 
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Since the orbital period and eclipse width are known, the masses and dimensions of the 
binary system can be computed once the mass ratio is specified (see Paper I). Unfortunately, 
the mass ratio of V Per has not been directly measured; however, its value can be estimated 
as follows. A lower limit on the mass ratio is established by noting that the observed eclipse 
width, A0 = 0.073 (Wood et al. 1992), requires that q (= M2/M1) ;^ 0.2 for orbital 
inclinations, i < 90° (Chanan et al. 1976). In addition, the requirement for stable mass 
transfer can be used to place an upper limit on the mass ratio given the mass of the secondary 
star. The dependence of M2 on q is weak, and to first order the mass of the secondary star 
can be expressed as a function of the orbital period alone. Prom Warner (1995, eqn. 2.100) 
we can write (M2/M0) ~ 0.065P^/^(/ir). For an orbital period of 2.57 hr, we can estimate 
M2 — 0.2Mq. Por secondary star masses less than 0.4 Mq, stable mass transfer requires 
g ^ 2/3 (e.g. Politano 1988). Thus, the mass ratio of V Per must lie roughly in the range 
0-2 ^ q ^ 0.6. Since for nova systems the mass of the white dwarf is almost certainly greater 
than 0.5 Mq, we can further reduce the upper limit of plausible mass ratios to g ~ 0.4, with 
the most likely value lying closer to g ~ 0.2, given that the average white dwarf masses in 
classical nova systems is beheved to be of order 0.85 Mq (Smith & Dhillon 1998). In order to 
explore the effect of the mass ratio on our analysis, we will explore solutions for mass ratios 
of both q — 0.2 and q — 0.4, with the former representing the most hkely model. PoUowing 
the analysis described in Paper I, the masses and dimensions of the V Per system have been 
computed for both of our representative mass ratios, and are summarized in Table 6. 

In addition to the mass ratio, the spectral type and temperature of the secondary star 
can be estimated prior to beginning the fitting procedure. In their statistical study of the 
properties of cataclysmic variable stars. Smith & Dhillon (1998) found an empirical relation 
between the orbital period and spectral type of the secondary star. Por an orbital period of 
2.57 hr, the most likely spectral type is between M4V and M5V. According to Popper (1980) 
a main sequence star with a spectral type in this range is expected to have a temperature 
of ~ 3300 K. Thus, for the purposes of our study, we will adopt T2 = 3300 K, and not vary 
this parameter during the fitting procedure. 

Another input parameter than can be constrained to a limited extent is the disk rim 
thickness, hj- = h/R^, where h is the height of the rim perpendicular to the plane of the 
disk and is the disk radius. According to the accretion disk models of Cannizzo (2001) 
the pressure scale height near the disk midplane gives hr — 0.01 to a good approximation. 
Unfortunately, this theoretical scale height, most often referred to by disk modelers, is dis- 
tinct from aspect ratios, /i/r ~ 0.03 — 0.1, that are measured by eclipse mapping techniques, 
and relevant to our study. In particular, Cannizzo & Wheeler (1984) carried out vertical 
integrations of accretion disk structure and showed that, even under standard a-disk theory, 
the disk photosphere can lie at ~ 3 — 4 pressure scale heights. According to Warner (1995, 
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eqn. 2.52), similar studies by Meyer & Meyer-Hofmeister (1982) and Smak (1992) show that 
the rim height can be approximated by h,- ~ O.OSSM^g^''. For plausible mass accretion rates 
between ~ 10^^ g s~^ and ~ 10^^ g s~^, we estimate rim heights of in the range of ~ 0.04 
to hj. ~ 0.08, which is consistent with our earlier estimate. To cover the range, we consider 
specific rim heights of both /ij. = 0.04 and /ir = 0.08 in our models of V Per. For comparisons 
we have also considered a = 0.0 "fiat disk" case for comparison with previous studies. The 
values of the remaining parameters, the "fitting parameters" , which include the radial disk 
temperature parameter a, and the temperatures of the disk rim, bright spot, and component 
stars, are constrained through the model fit. 



5. Light-Curve Fitting 

To begin the fitting procedure, we select plausible ranges of values for the five fitting 
parameters, indicated in Table 5 as "variable" . Then each parameter is varied one at a time 
through a suitable range while the remaining parameters are held constant. In this way a 
5-dimensional parameter space of models is explored for plausible solutions. The goodness 
of fit at each stage is determined using a standard ^Gst. The deviations between the 
model and the data are determined for each of the four colors, and the final statistic is 
determined by weighting each of the colors equally. Since we model the light curve near 
eclipse only, the calculation of has been restricted to orbital phases between = —0.12 
and (p = 0.12, which corresponds roughly to the time just prior to the onset of disk ingress to 
the completion of bright-spot egress. The range over which each parameter is varied has been 
chosen to be sufficiently large so that all plausible model solutions were explored. A thorough 
sampling required eight values to be considered for each fitting parameter. Thus, a total of 
32,768 (= 8^) models were computed for each mass ratio and rim thickness considered. The 
best-fitting solutions for our grid of models are given Table 7. 

In their earlier monochromatic ( "white light" ) eclipse-mapping study of V Per, Wood 
et al. (1992) found that the brightness temperature gradient of the accretion disk could be 
made consistent with a steady-state {a — 0.75) temperature profile only if the inner edge 
of the accretion disk was of order it!in ~ 0.2 i?Lij where i?Li is the distance from the center 
of the white dwarf to the inner Lagrangian point. These authors pointed out that such a 
situation might arise if the inner accretion disk was disrupted by the white dwarf's magnetic 
field, as is believed to be the case in intermediate polar systems (e.g. Warner 1995). To 
explore this possibility further, we have computed an additional set of eclipse models where 
the inner disk radius has been set at the value, i?in = 0.2i?Li, adopted by Wood et al. (1992). 
The results of our intermediate polar models are shown in parentheses in Table 7. 



5.1. The Model Solution Distributions 



Exploring the matrix of possible solutions reveals that many combinations of parameters 
produce good fits to the data. The goal is to determine which particular combination of 
parameters best represents the true model for V Per. The best model fits are characterized 
by a reduced chi-square, xt — 1-5; however, a wide range of parameters will give comparable 
fits to the light curves. To explore the range of plausible solutions, frequency distributions for 
the top 100 best-fitting model solutions were constructed for each of the fitting parameters. 
Figures 3-5 show the frequency distributions of the five fitting parameters, together with 
the the model B — V color, for our preferred mass ratio oi q = 0.20. Figure 3 shows 
the distributions for our rimless model with = 0.0, with Figures 4 and 5 showing the 
distributions for /ij. = 0.04, and /ir = 0.08, respectively. The mean and standard deviation 
of each parameter distribution, along with the mean xl ioi each model, are summarized in 
Tables 8 and 9 for the standard and truncated inner disk models, respectively. The mean 
xl of the top 100 model solutions is significantly higher in the rimless disk models, but still 
refiects acceptable fits to the data. As long as the number of model solutions included in 
the distributions is sufficiently small so that <xl> ^ 1-8 (which represents a plausible fit to 
the data), the distribution means do not depend strongly on the precise number of model 
solutions included in the distributions. 

The frequency distributions reveal how well a given parameter is constrained by the 
model fit. If a parameter is tightly constrained by the data, its frequency distribution will 
be narrow and centered on its optimum vahie. Thus, the mean and standard deviation of 
each distribution not only provide an estimate of the optimum value for a specific param- 
eter, but also a quantitative assessment of its uncertainty. The most tightly constrained 
parameter is clearly the disk radius, R^, where a radius of ~ 70 — 75% of the distance to 
the inner Lagrangian point is required by virtually all the models. The value of this pa- 
rameter is primarily determined by the phase width where the disk eclipse begins and ends, 
and is relatively insensitive to other system parameters (e.g. Sulkanen et al. 1981). The 
light curves do not exhibit any significant asymmetry of the eclipse profile, and there is no 
prominent "hump" prior to eclipse, so it is not surprising that none of the models considered 
required significant radiation from the bright spot in the system. Furthermore, the eclipse 
profile, which lacks the sharp transitions on ingress and egress characteristic of the white 
dwarf eclipse suggest that the white dwarf star does not contribute significantly to the sys- 
tem light. Thus, as expected, none of the models we considered required a luminous (hot) 
white dwarf. As an example. Figure 6 shows the best-fitting q = 0.20 and hr — 0.08 model 
plotted together with the phased echpse data. The best fitting models for q — 0.4, and for 
other rim thicknesses fit the data almost as well, and cannot be ruled out based on the fit 
alone. 



- 11 - 



6. Discussion 

6.1. Disk Structure 

Monochromatic eclipses of V Per have been studied by Wood et al. (1992) using a 
maximum entropy echpse-mapping technique. A principal conclusion of their study was 
that the brightness temperature distribution of the accretion disk was flatter than that 

expected for a steady-state disk. These authors showed, however, that if the inner disk in 
V Per were disrupted out to ~ 0.15 — 0.25i?Lii as might be expected if V Per harbors 
a magnetic white dwarf, then a brightness temperature gradient very close to that expected 
for a steady-state disk could be achieved. 

Our multi-color eclipse simulations are qualitatively in agreement with the results of 
Wood et al. (1992). Referring to Table 7, we see that for a given mass ratio and rim 
thickness, the intermediate polar models are generally characterized by slightly larger values 
of the disk temperature parameter, a. However, in our case, only models with a disk rim 
produced a best-fitting value of a close to the steady-state value. Only two cases - the q — 0.2, 
hr — 0.08 standard model and the q — 0.4, — 0.08 intermediate polar model - resulted 
in a value of a = 0.75. It is noteworthy that for a given rim height the value of a is 
systematically higher in all of the intermediate polar models compared with the standard 
models, underscoring the synergistic effect between the rim and the truncated inner disk in 
achieving a steady-state disk solution. For the orbital inclination of 85.4° required by our 
g = 0.2 model, the central region of the accretion disk, including the white dwarf, is just 
barely visible with a rim thickness of h,- = 0.08. Any "thick disk" models characterized by 
rim height > 0.08 would result in complete shadowing of the white dwarf and part of the 
disk at all orbital phases, and thus would be insensitive to the white dwarf's luminosity. 

Wood et al. also considered an alternative model with a disk rim, and found that their 
rim model did not result in a significantly steeper brightness temperature gradient. However, 
the rim model they considered was considerably less sophisticated than ours. In particular, 
their models employed a "dark" disk rim solely to study the effect produced by shadowing 
of the disk by the rim. The effect of the rim's contribution to the overall disk hght was not 
included in their models. Furthermore, rim heights of only 1-2% of the disk radius were 
considered, which required inclination angles, i > 89°, to achieve significant rim shadowing. 

As a test of the importance of a bright rim versus a dark rim like that considered by 
Wood et al., we ran a series of models with our standard flared disk, but with the rim 
brightness switched off. These "dark rim" models generally produced the same best-fit 
solutions as our standard models, but with one difference: a significantly lower value of 
the disk radial temperature parameter, a. Specifically, for our four standard models with 
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a rim {q = 0.2, = 0.04; g = 0.2, K = 0.08; g = 0.4, /i, = 0.04; g = 0.4, /i, = 0.08), 
we found values of a = 0.15, 0.45, 0.25, and 0.35, respectively, compared with the values 
of a = 0.35, 0.75, 0.25, and 0.45 from our bright rim models (see Table 7). The only case 
where the bright rim apparently made no difference was in the model characterized by 
q = 0.4, hj: = 0.04. However, in this case, the lower inclination and modest rim height make 
the contribution of the rim to the total disk light the smallest of all models considered. 

The importance of a bright disk rim can be best visualized by referring to Figures 7 
and 8, which show the system geometry at mid-eclipse. At high inclinations the disk rim 
makes a significant contribution to the total disk light. Since the disk rim (excluding the 
bright spot) has a nearly constant surface brightness, the eclipse of this rimmed disk will be 
more V-shaped than the eclipse of a flat disk, even in the case where the disk has a steady- 
state {a — 0.75) radial temperature profile (e.g., see the eclipse simulations and discussion of 
DW UMa by Knigge et al. 2000, 2004). As mentioned earlier, Smak (1994) has shown that 
modeling such eclipses under the assumption that the disk is fiat (as in eclipse mapping) 
leads to the spurious conclusion that the disk temperature profile is shallower that a = 0.75. 
In the case of V Per, an additional effect is that for both mass ratios considered, the accretion 
disk is never completely eclipsed. In models with a disk rim, the contribution of the rim 
to the total system light at mid-eclipse can be significant. In order to produce an eclipse 
with the observed depth (in the normalized light curves), the inclusion of a bright outer 
disk rim requires that the eclipsed inner disk be sufficiently bright to compensate, i.e. to 
keep the ratio of the mid-eclipse to out-of-eclipse system brightness unchanged. For a given 
outer disk temperature, this can only be accomplished by having a luminous white dwarf, 
or by having a steep disk temperature profile, or both. These two possibilities produce very 
different ingress and egress eclipse profiles, and can therefore be distinguished by the model 
fit. In the case of V Per, the steep ingress and egress characteristic of a significant white 
dwarf contribution to the overall luminosity is not observed, leaving a steep disk temperature 
gradient as the only viable alternative. 



6.2. Correlations Between Model Parameters 

Further insights into the the effect that the various input parameters have on acceptable 
model solutions can be found by exploring correlations between the parameters. We begin by 
considering an expanded set of model solutions defined by xl < 2.0, which allows potential 
correlations to be studied over an expanded range in parameter space. Correlations between 
pairs of parameters are explored by allowing these parameters to vary while holding the 
remaining parameters fixed at their optimum values. For a total of five fitting parameters. 
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there are a total of 10 such pairings. As a representative example, the correlations for the 
best-fitting q — 0.2, h^. — 0.04 model are shown in Figure 9. 

As in previous studies of GY Cnc and EX Dra (Papers I and II), the most highly 
correlated pairs involve i?d, T^, and a, which together determine the disk luminosity. There 
is a major difference in the correlations found in the present study, however. In the cases of 
GY Cnc and EX Dra, which are longer orbital period systems with larger mass ratios, the 
accretion disks were entirely occulted at mid-eclipse. Thus, the depth of eclipse fixes the disk 
luminosity. For a constant disk luminosity, the radial temperature parameter a should be 
negatively correlated with the temperature at the disk periphery, Td, and the disk radius, /?d, 
as was indeed observed for these two systems. In the case of V Per, however, the situation is 
more complicated as described in the previous section. Because the disk is never completely 
obscured, the eclipse depth no longer uniquely defines the disk luminosity. The normalized 
eclipse depth can be achieved with a relatively low- luminosity disk having a relatively shallow 
radial temperature gradient, or with a disk having a bright rim (relatively high temperature 
and/or vertical extent) and steep temperature gradient (large a). As the contribution of 
the disk rim is increased, the inner disk luminosity must increase to compensate, requiring 
a steeper temperature gradient and a higher value of a. Thus, for a given disk radius and 
rim thickness, a should be positively correlated with the outer disk (rim) temperature, T^, 
as observed in the case of V Per, and shown in Figure 9. In view of these considerations, it 
is not surprising that the q — 0.40, — 0.00 and q — 0.20, — 0.08 models yield the least 
luminous and most luminous disks, respectively. 

In principle, a smaller disk with the appropriate luminosity that is completely eclipsed 
could produce the observed eclipse depth, however, the shape of the eclipse profile, not just 
its depth is important for a good model fit. In the case of V Per, the maximum width of 
the eclipse profile (as defined by the onset of ingress and the termination of egress) requires 
a relatively large disk radius that approaches the tidal limit (see Figure 6) . 

6.3. Model Constraints from the Disk Luminosity 

Although a range of models can produce acceptable fits to the data, as noted above, 
these models differ in one important respect: they represent significantly different disk lu- 
minosities. Referring to our best- fitting model solutions (Table 7), we are reminded that 
rimless models are characterized by flat temperature profiles (small a). For a given outer 
disk radius and temperature, disks with relatively flat temperature gradients will be con- 
siderably less luminous than steady-state (rimmed) disks characterized by a = 0.75. Thus, 
a potential way to constrain the models further is to favor those models that produce the 
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expected disk luminosity, and absolute magnitude for V Per. The disk luminosities can be 
explored through two independent approaches: directly from our model fluxes, and from an 
estimate of the mass accretion rate inferred from our models. 

In the first approach, we can compare the ratio of our model V fluxes for the disk and 
secondary star as viewed from an inclination angle of i — (i.e. "face-on"). This quantity, 
7/d.s = -^v(^ = ^)/Fy{i — 0), is given in Table 10 for our best-fitting models. If we ignore 
limb-darkening and assume a spherical secondary star and a flat circular disk, then the 
luminosity ratio is approximately equal to half the model flux ratio (a flat disk has half the 
surface area of a sphere). Thus, Ly/Ly ~ ?7d,s/2, and an estimate of the secondary star's 
luminosity will yield a crude estimate for the disk luminosity. In reality, this estimate will be 
a slight underestimate of the disk luminosity because we have neglected the relatively small 
contribution of the rim to the overall disk luminosity. As an example, the flux ratio of our 
best-fitting q — 0.2, — 0.08 model is rjd^s — 3.18 x 10^, which represents a luminosity ratio 
of ~ 1.59 X 10^. Thus, we estimate the disk in this model to be ~ 8.0 mag brighter than the 
secondary star. An estimate of My (2) can be found given the secondary star's spectral type, 
which wc estimated earlier to be ~M4V-M5V. For a main-sequence secondary star of radius 
7?2 = 0.24 Rq (Table 6) and spectral type M4V-M5V, we find 12.4 < My (2) < 12.8 (Popper 
1980). Thus, for this model, we estimate an absolute visual magnitude for the accretion disk, 
and in effect for the V Per system (the white dwarf and secondary star contribute little V 
hght compared to the disk), to be My ~ 4.6. Although somewhat bright for a P = 2.6 hr 
system, this value is certainly plausible for V Per, and is typical for the absolute visual 
magnitude of nova remnants (e.g. Warner 1995). 

As a check on this result, we can also make a rough estimate the disk luminosity from 
an estimate of the mass transfer rate coupled with the M vs. My relation from Tylenda 
(1981). We begin with a well-known expression for the effective temperature of the accretion 
disk (e.g. Pringle 1981): 



where it!in is the inner radius of the disk. By adopting the model radius and temperature of 
the outer disk (rim) , we can compute the mass accretion rate for a given model. Specifically, 




(10) 




1 



gm s 



1 



(11) 



Mass transfer rates (M17 = M/10^^ g s~^) for our best-fitting models are given in Table 10. 
For the specific case of the q — 0.2, hr — 0.08 model with a — 0.75, the mass-transfer rate 
can be converted to an equivalent My through an extrapolation of the relation given in 
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Table 1 of Tylenda (1981).^ We find a value of My ~ 4.3, which is consistent with the value 
of My ^ 4.6 estimated earher from the disk-to-secondary star luminosity ratio. 

The absolute magnitudes and accretion rates of systems with fiat accretion disks {My ~ 
7 — 8 and M^ ~ 1.5 — 2.5) are more typical of dwarf nova disks at minimum light. Such 
feeble disks with their shallow temperature gradients {a — 0.15) are unlikely to be seen in 
nova and nova-like systems with stable accretion unless a significant amount of mass and 
accretion energy is being lost from the inner disk, for example, by driving a wind. In such a 
case, the mass accretion rate in the outer disk could possibly be increased without necessarily 
producing a luminous disk with a steep temperature profile. 

Since V Per does not exhibit dwarf nova eruptions, a further constrain on possible 
models for V Per is provided by the requirement that the mass transfer rate exceed the 
citical value for stable accretion. Merit- A convenient expression for the critical mass transfer 
rate is given by Merit — 10^^ Ril^^ M^ gm s~^, where Rio is the radius of the disk in units 
of 10^° cm (e.g. see Shafter et al. 1986, and references therein). Values of Merit for our 
best-fitting models are given in Table 10. The fiat disk models have mass accretion rates 
that are significantly below the critical value, while only the thick disk models have rates 
clearly in excess of the critical value. 

To minimize the number of fitting parameters in our models, we have set the disk 
rim temperature equal to the minimum temperature in the outer disk, T^. It seems possi- 
ble, however, that this assumption overestimates the value of used to compute the disk 
temperature profile since energy input from the mass-transfer stream is likely to raise the 
temperature of the disk rim over that appropriate for the outer disk region. If so, then our es- 
timates of mass transfer rates and disk luminosities are best considered upper limits. In this 
case, the flat disk models with shallow temperature gradients become even more problematic, 
and our conclusion favoring the high- mass-accretion-rate, high-luminosity, thick-disk model 
for the V Per system is strengthened. 

6.4. The Distance to V Per 

The distance to V Per can be estimated by comparing the observed brightness at mid- 
eclipse with estimates of the absolute magnitude of the secondary star, and the fraction of 



^It has been necessary to apply a small correction (~ 0.8 mag) to the resulting absolute magnitude to 
account for the difference in Mi and Rd between Tylenda's model (Mi = IMq, = 5 x 10^° cm) and our 
model. Following Warner (1995), the correction for R^ has been estimated from the models of Wade (1984). 
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light it contributes at mid-eclipse. The fraction of light contributed by the secondary star at 
mid-eclipse, /2 = /v(0 = O)//v'*(0 = 0), depends on the model considered, and in particular 
on the properties of the disk rim. Values of /2 for each model are given in Table 10. In the 
rimless models, the disk contributes relatively little to the light at mid-eclipse compared with 
the models with rims. According to our best-fitting solutions for models with a thick rim, the 
secondary star contributes ~ 1.0% of the light at mid-eclipse, while in the rimless models, 
where the disk luminosity is relatively low, the secondary star contributes a substantial 
~ 20% of the mid-echpse light. Given that my — 19.34 (see Table 3) at mid-echpse, we 
estimate the apparent distance moduli (m — M)v given in Table 10. 

The interstellar extinction along the line of sight can be estimated from the color excess 
obtained from a comparison of the observed and model B — V colors. V Per has an observed 
color of S — y ~ 0.6 ± 0.2 outside of echpse, where much of the uncertainty is due to 
the limited phase coverage of our outside echpse S-band data. Assuming a ratio of total-to- 
selcctivc extinction of Ry = 3 yields the visual extinctions. Ay, along the line of sight to V per 
given in Table 10. Our models with a rim have colors outside of eclipse that are significantly 
bluer than the colors of the rimless models, leading to larger values of Ay. Correcting the 
apparent distance moduli yields the model-dependent distance estimates given in Table 10. 

As a consistency check, the distance to V Per can also be determined from a comparison 
of the out-of-eclipse magnitude with an estimate of the system's absolute magnitude corrected 
for inclination. Warner (1987, 1995) gives a correction factor for the conversion of absolute 
magnitude to apparent absolute magnitude of AMy{i) = —2.5 log [(1 + 1.5 cos i) cos i]; 
however, this correction does not take into account any contribution from the disk rim.^ 
Fortunately, the ratio of rim-to-disk light (r^r.d) is given by our models (see Table 10). In the 
case of our favored q — 0.2, hr — 0.08 model, the rim contributes ~ 50% of the disk light. 
Applying these corrections for inclination, taking My — 4.6 for the q — 0.2,h,- — 0.08 model, 
and adopting my — 18.1 for V Per out of echpse (see Table 3), yields a reddening-corrected 
distance modulus of (m — M)y = 9.6, and a corresponding distance of ~ 820 pc. Given the 
uncertainties, this value is consistent with our estimate of 1 kpc based on the properties of 
the secondary star. 

In summary, we estimate that V Per lies at a distance of ~ 1.0 kpc, in our favored 
scenario where the disk has a significant rim that contributes to the light at mid-eclipse. 
In contrast, our rimless models, which predict that a significantly larger fraction of the 
light at mid-eclipse arises from the secondary star, suggest that V Per lies at a distance of 



•^The expression given in Warner (1995) corrects an error in the Warner (1987) formula, however the sign 
is incorrect for his stated conversion of apparent absolute magnitude to absolute magnitude. 



-17- 



only ~ 0.5 kpc. In both cases, a significant component of the uncertainty in our distance 
estimates arises from uncertainties in the visual extinction. Overall, our favored distance 
estimate is consistent with the value of ~ 1 kpc derived by Shatter (1997) from a comparison 
of the apparent magnitude at maximum light (mmax — 4), an assumed absolute magnitude 
at maximum light. My = —7.5, and a visual extinction of Ay = 1.4 estimated from the 
reddening maps of Neckel et al. (1980). 

7. Conclusions 

We have performed the first multicolor eclipse study of the classical nova V Per using 
a parameter-fitting eclipse code to model the light curve. Our study has explored a variety 
of assumptions regarding the structure of the accretion disk, including models with a large 
inner disk radius, as might be expected if V Per has a strongly magnetic white dwarf, and 
models that include a disk rim. In agreement with the earlier monochromatic maximum 
entropy eclipse mapping study of V Per by Wood et al. (1992), we find that a fiat accretion 
disk that extends down to the surface of the white dwarf must have a brightness temperature 
profile that is significantly fiatter than that expected for a steady-state disk. In an attempt 
to resolve the discrepancy. Wood et al. explored models where the inner disk was disrupted, 
as might be expected if V Per harbors a magnetic white dwarf. They found that temperature 
profiles consistent with the steady-state case could be produced in models with an increased 
inner disk radius, i?in ~ 0.2i?Li- 

Although the brightness temperature profiles of our best-fitting models become steeper 
when the inner disk is truncated out to Ri^ = 0.2i?Lij unlike Wood et al., we were imable to 
reproduce a flat-disk model with a ~ 0.75 with a truncated inner disk alone. However, our 
models that included a disk rim were characterized by considerably steeper brightness tem- 
perature profiles, particularly in our "thick disk" (hr — 0.08) case where we found a — 0.75 
gives the best fit to the data. Although acceptable fits for models with rims can be found 
with smaller values of a and lower outer disk temperatures (faint rims), the only way to 
achieve a — 0.75 is with a model having a bright rim. Furthermore, only the bright rim 
models give accretion rates and absolute magnitudes that are typical of nova systems with 
stable accretion. By extrapolating our result for V Per, we confirm the conclusion reached 
by Smak (1994) from his generic eclipse simulations, that prior attempts to analyze cata- 
clysmic variable eclipses with flat-accrction-disk models have significantly underestimated 
the brightness temperature gradients of the accretion disks. 

In summary, regardless of the adopted mass ratio or rim thickness, we find that V Per 
is most likely characterized by a large accretion disk (extending to near its tidal radius). 
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characterized by Rd c::^ 0.75i?Li, with a temperature Td ~ 6000 — lOOOOK at its outer edge. 
No evidence was found for an unusually luminous (hot) white dwarf or bright spot in the 
system. However, if the disk is in a steady-state, optically-thick configuration, we have 
shown that the disk must possess a significant rim. As in V Per, it is also likely that other 
short-period cataclysmic variable stars with mass ratios, q <^ 0.5, will also potentially harbor 
accretion disks that are not fully occulted at mid-eclipse. As shown in the present study, it 
is imperative that any attempts to model the accretion disks in these systems should take 
into account the effect of a disk rim. 

We thank the staff of the Steward Observatory for granting time on the Bok 2.3 m 
reflector and Karl Gordon for assistance with the observations. We also thank John Cannizzo 
and Ted Daub for reading an earlier version of the manuscript and for subsequent discussion, 
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Table 1. CVs in the Period Gap 



Name 


lype 


Orbital Period (hr) 


V516 Pup 


AM 


2.29 


SDSS J205017-053627 


NL 


2.30 


QS Tel 


AM 


2.33 


DD Cir 


N, DQ 


2.34 


NY Ser 


SU 


2.35 


SDSS J080908+381406 


NL 


2.38 


V348 Pup 


DQ 


2.44 


MN Dra 


SU 


2.50 


RX J1554.2+2721 


AM 


2.53 


V Per 


N 


2.57 


V795 Her 


DQ 


2.60 


IRXS J052430+424449 


AM, DQ 


2.62 


V349 Pav 


AM 


2.66 


QU Vul 


N 


2.68 


SDSS J075240+362823 


AM? 


2.70 
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Table 2. Summary of Observations 





UT Time 


Time Resolution* 


Number of 




UT Date 


(start of observations) 


(sec) 


Exposures 


Filter 


2002 Oct 10 


4:57:00 


41.5 


371 


V 


2002 Oct 10 


10:34:00 


41.5 


151 


V 


2002 Oct 11 


7:05:00 


57.0 


248 


B 


2002 Oct 29 


4:33:00 


56.1 


74 


R 


2002 Oct 29 


7:00:00 


56.1 


74 


R 


2002 Oct 29 


9:11:38 


56.2 


121 


I 


2002 Dec 09 


2:20:35 


56.4 


111 


I 


2002 Dec 09 


5:22:00 


56.4 


77 


B 


2002 Dec 09 


7:32:28 


56.4 


96 


I 


2002 Dec 11 


5:17:00 


56.4 


132 


R 



^Mean time interval between exposures (Integration time plus readout time) 
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Table 3. Mean Magnitudes and Colors 



Effective 
Photometric Wavelength 



Parameter (A) Out of Eclipse*^ Primary Minimum 



B 4375 18.7 ±0.1 20.05 ± 0.09 

V 5460 18.1 ±0.1 19.34 ±0.08 

R 6543 18.1 ±0.2 19.13 ±0.31 

1 7992 17.5 ±0.1 18.61 ±0.17 

B-V... ... 0.6 ±0.2 0.71 ±0.12 

V-R... ... 0.0 ±0.3 0.21 ±0.32 

V-I... ... 0.6 ±0.2 0.73 ±0.19 



^The mean magnitude of the hght curve excluding orbital phases in 
the range: -0.15 < < 0.15. 
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Table 4. Eclipse Timings 



HJD (mid-cclipse) 


Cycle Number 




0-C 




(2,400,000+) 


(E) 


Filter 


(xlO^ day) 


reference 


47442.8264. . . 


0. 


BVR 


0.438541 


1 


47442.9337. . . 


1. 


BVR 


0.615067 


1 


47443.7896. . . 


9. 


BVR 


-0.472726 


1 


47443.8965. . . 


10. 


BVR 


-0.696200 


1 


47444.0040. . . 


11. 


BVR 


-0.319674 


1 


47446.8969. . . 


38. 


BVR 


0.246525 


1 


47447.0038. . . 


39. 


BVR 


0.023051 


1 


47448.9321. . . 


57. 


BVR 


0.100518 


1 


50430.7140. . . 


27892. 


V 


0.099267 


2 


50430.8207. . . 


27893. 


V 


-0.324207 


2 


50481.5978... 


28367. 


R 


0.249075 


2 


50481.7052... 


28368. 


B 


0.525601 


2 


50509.6639. . . 


28629. 


V 


-0.001136 


2 


50692.8455. . . 


30339. 


B 


0.458173 


2 


50695.8445. . . 


30367. 


R 


0.000898 


2 


50696.9157. . . 


30377. 


B 


-0.033843 


2 


50696.8089. . . 


30376. 


R 


0.289631 


2 


50730.7668. . . 


30693. 


I 


0.048345 


2 


50730.8740. . . 


30694. 


I 


0.124871 


2 


50730.9807. . . 


30695. 


I 


-0.298603 


2 


50731.9451. . . 


30704. 


V 


-0.009870 


2 



Table 4 (cont.). Eclipse Timings 



HJD (mid-eclipse) Cycle Number O — C 



(2,400,000+) 


{E) 


Filter 


(xlO^ day) 


ref 


oO/o/.obyi. . . 


or\c\ACi 
o0y4D. 




0.109401 


o 
z 


oO/oo. /zoo. . . 


or\n A 
o0yo4. 


T 
1 


— O.D/Ooyz 


o 
z 


50758.8331... 


30955. 


B 


-0.001866 


2 


50759.7971... 


30964. 


B 


-0.113133 


2 


50759.9041... 


30965. 


V 


-0.236607 


2 


52172.541... 


44152. 


R 


-0.589411 


3 


52557.7580. . . 


47748. 


V 


0.397768 


4 


52557.8642. . . 


47749. 


V 


-0.525707 


4 


52557.9720. . . 


47750. 


V 


0.150819 


4 


52558.8293. . . 


47758. 


B 


0.463027 


4 


52558.9362. . . 


47759. 


B 


0.239553 


4 


52576.7186. . . 


47925. 


R 


0.142854 


4 


52576.8255. . . 


47926. 


R 


-0.080620 


4 


52576.9323. . . 


47927. 


I 


-0.404094 


4 


52617.6398... 


48307. 


I 


0.175752 


4 


52617.7469. . . 


48308. 


B 


0.152278 


4 


52617.8537. . . 


48309. 


I 


-0.171196 


4 


52619.7820. . . 


48327. 


R 


-0.093730 


4 



^1.- Shatter & Abbott 1989; 2.- Dahm (1997); 3.- Katysheva et al. 
2002; 4.- This work. 
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Table 5. Model Input Parameters 



Parameter Definition Value^ 



P Orbital Period, P 2.57 hr 

q Mass ratio, M2/M1 0.20,0.40 

A(f) Eclipse phase width 0.073 

i Orbital inclination computed 

a Disk temperature parameter variable 

i?d Disk radius variable 

i?in Inner disk radius Ri, 0.2i?Li 

hr Disk rim parameter 0.00,0.04,0.08 

Td Temperature of disk perimeter variable 

Ri Radius of white dwarf computed 

i?2 Radius of secondary star computed 

Ti Temperature of white dwarf variable 

T2 Temperature of secondary star 3300 K 

Xs Bright spot temperature factor variable 

Rg Bright spot radius 0.2it!d 



^Where given explicitly, the values are fixed and not varied 
during the fitting procedure. The values of "variable" parame- 
ters are determined by the model during the fitting procedure, 
while the "computed" values are calculated by the model for a 
given mass ratio. 
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Table 6. Binary Parameters 



Parameter 


q = 0.20 


q = 0.40 


in 


85.4 


79.4 


a{Re) 


0.95 


0.80 


Rhiia) 


0.66 


0.59 


i?i(10-2 Rq)... 


0.97 


1.52 


M,{Mq) 


0.85 


0.43 


R2{Rq) 


0.24 


0.24 


M2{Mq) 


0.17 


0.17 



Table 7. Model Grid Solutions^ 









/ 


0.20 










0.40 






Parameter 


K = 0.00 


hr — 


: 0.04 


hr = 0.08 


hr — 


: 0.00 


hr — 


: 0.04 


hr - 


= 0.08 


Rd{a) 


0.49 (0.40) 


0.49 


(0.40) 


0.46 (0.40) 


0.41 


(0.41) 


0.44 


(0.44) 


0.41 


(0.38) 


-Rd(a)lim'^ 


0.50 (0.50) 


0.50 


(0.50) 


0.50 (0.50) 


0.43 


(0.43) 


0.43 


(0.43) 


0.43 


(0.43) 


RdiRLi) 


0.65 (0.60) 


0.75 


(0.65) 


0.75 (0.70) 


0.70 


(0.65) 


0.75 


(0.70) 


0.75 


(0.75) 


a 


0.15 (0.35) 


0.35 


(0.55) 


0.75 (0.75) 


0.15 


(0.25) 


0.25 


(0.45) 


0.45 


(0.65) 


71(10^ K) 


10.0 (10.0 ) 


10.0 


(15.0) 


10.0 (25.0) 


10.0 


(10.0) 


10.0 


(10.0) 


10.0 


(15.0) 


7d(103 K) 


6.0 (6.0) 


6.0 


(6.0) 


10.0 (6.0) 


6.0 


(6.0) 


6.0 


(6.0) 


8.0 


(6.0) 


Xs 


1.2 (1.2) 


1.2 


(1.1) 


1.3 (1.1) 


1.1 


(1.1) 


1.3 


(1.2) 


1.3 


(1.2) 


{B - V). 


0.63 (0.56) 


0.45 


(0.44) 


0.03 (0.34) 


0.59 


(0.56) 


0.49 


(0.45) 


0.19 


(0.33) 


Xi/,min 


1.7 (1.7) 


1.6 


(1.6) 


1.5 (1.6) 


1.6 


(1.6) 


1.5 


(1.6) 


1.5 


(1.5) 



'^Values in parenthesis represent intermediate polar models where the inner accretion disk has 
been truncated at RdiRhi) = 0.2. 

^The limiting tidal radius of the disk computed using eqn (2.61) of Warner (1995). 



Table 8. Standard Model Distribution Mean Solutions 







q = 0.20 






q = 0.40 




Parameter 


K = 0.00 


K = 0.04 


K = 0.08 


K = 0.00 


K = 0.04 


K = 0.08 


<xl> 


1.74 


1.61 


1.59 


1.74 


1.66 


1.58 




0.67 ±0.05 


0.73 ±0.04 


0.73 ±0.03 


0.71 ±0.04 


0.73 ±0.03 


0.75 ±0.02 


a 


0.17 ±0.09 


0.45 ±0.13 


0.60 ±0.13 


0.19 ±0.10 


0.32 ±0.11 


0.51 ±0.12 


Ti{10^ K) 


14.8 ±2.05 


23.0 ± 4.44 


25.6 ±4.77 


14.7 ±4.52 


18.5 ±4.51 


24.0 ±5.84 


Td(lO^ K) 


7.22 ± 0.97 


8.76 ± 1.91 


8.92 ± 1.77 


8.15 ± 1.68 


8.35 ± 1.82 


10.2 ±2.40 


Xs 


1.18 ±0.15 


1.27±0.11 


1.23 ±0.08 


1.14 ±0.14 


1.29 ±0.13 


1.32 ±0.11 




0.35 ±0.15 


0.13 ±0.19 


0.08 ±0.14 


0.26 ±0.17 


0.18 ±0.19 


0.03 ±0.16 



Table 9. Intermediate Polar Distribution Mean Solutions 







q — 0.20 






q — 0.40 




Parameter 


K = 0.00 


hr = 0.04 


hr = 0.08 


hr: = 0.00 


hr: = 0.04 


hr = 0.08 


<xl> 


1.77 


1.67 


1.69 


1.73 


1.65 


1.62 




0.65 ±0.05 


0.68 ±0.05 


0.65 ±0.05 


0.69 ±0.05 


0.72 ±0.04 


0.72 ±0.03 


a 


0.42 ±0.15 


0.60 ±0.14 


0.60 ±0.15 


0.35 ±0.15 


0.55 ±0.16 


0.65 ±0.12 


Ti{10^ K) 


15.2 ±2.33 


23.8 ±4.11 


24.7 ±4.44 


13.2 ±2.57 


17.8 ± 3.02 


24.1 ± 4.42 


Td(lO^ K) 


6.54 ±0.88 


6.66 ± 0.94 


6.53 ±0.88 


6.82 ± 0.98 


6.94 ± 1.10 


7.36 ± 1.31 


Xs 


1.18 ±0.15 


1.16 ±0.08 


1.10 ±0.06 


1.15 ±0.14 


1.20 ±0.13 


1.18 ±0.10 




0.30 ±0.10 


0.26 ±0.13 


0.30 ±0.12 


0.31 ±0.14 


0.23 ±0.16 


0.18 ±0.14 



Table 10. Luminosities, Mass Accretion Rates, and Distances 







q = 0.20 








Q = 


0.40 




Parameter 


hr = 0.00 


hr = 0.04 


h, = 0.08 




0.00 


hr — 


: 0.04 


hr = 0.08 




155 (148) 


387 (283) 


3178 (493) 


98.4 


(83.4) 


148 


(146) 


622 (277) 


'/r,a 


fO 0) 


34 (0 34) 


51 (0 40) 


0.0 


(0.0) 


0.21 


(0 17) 


31 (0 19) 


My 


7.9 (7.9) 


6.9 (7.2) 


4.6 (6.6) 


8.4 


(8.6) 


7.9 


(7.9) 


6.4 (7.2) 


Mir (g s-i) 


1.49 (2.34) 


2.26 (2.83) 


17.4 (3.38) 


1.68 


(2.37) 


2.05 


(2.84) 


6.47 (3.36) 


Mcrit,17 (g S-^) 


1.72 (2.34) 


2.50 (1.72) 


2.50 (2.09) 


1.82 


(1.50) 


2.20 


(1.82) 


2.20 (2.20) 


/2 


0.23 (0.26) 


0.077 (0.11) 


0.0089 (0.056) 


0.16 


(0.19) 


0.088 


(0.096) 


0.022 (0.051) 


(m - M)v 


8.34 (8.20) 


9.52 (9.14) 


11.87 (9.87) 


8.73 


(8.54) 


9.38 


(9.29) 


10.88 (9.97) 


Ay 


0.0 (0.12) 


0.45 (0.30) 


1.71 (0.78) 


0.03 


(0.12) 


0.33 


(0.45) 


1.23 (0.81) 


d (kpc) 


0.47 (0.42) 


0.65 (0.59) 


1.08 (0.66) 


0.55 


(0.48) 


0.65 


(0.59) 


0.85 (0.68) 
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Fig. 1. — The multicolor light curves of V Per. For clarity of presentation, the data have 
been offset by addition of the following constants: i^-band: C=0.0; l^-band: C=1.2; i?-band: 
C=2.4; /-band: C=4.0. 
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Fig. 2. — The residuals of the the observed times of mid-echpse with respect to eqn (1) is 
plotted as a function of cycle number. The final grouping of points near cycle number 48000 
are based on eclipse timings from this work, while the remaining points are based on timings 
from the literature. There is no evidence for any period change in V Per over the ~ 15 years 
of available eclipse timings. 
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Fig. 3. — The q — 0.20, — 0.00 frequency distributions for each fitting parameter. The 
range of each parameter has been chosen to include all combinations of parameters that 
result in acceptable fits of the top 100 best-fitting models to the data. The cross-hatched 
regions indicate the parameter values that produce the optimum fit to the data. 
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Fig. 4. — The same as Figure 3, but for h^. = 0.04. 
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Fig. 5. — The same as Figure 3, but for = 0.08. 
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Fig. 6. — The best-fitting eclipse profiles (solid lines) for the q = 0.20, /ir = 0.08 model 
are plotted together with the observed data. The data have been converted to relative 
intensity and normalized to unity in the orbital phase range: 0.15 < (f) < 0.20. For clarity 
of presentation the V, R, and /-band light curves and models have been shifted upward by 
constant offsets of 1, 2, and 3, respectively. 
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Fig. 7. — A schematic representation of the V Per geometry at mid-echpse for the q = 
0.2, i = 85.4, /ij. = 0.04 modeL Note that the accretion disk is never completely eclipsed, and 
that the outer disk and rim make a significant contribution to the light at mid-eclipse. 
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Fig. 9. — The correlations between pairs of model parameters for the representative q = 
0.20, /ir = 0.04 model are shown for each of the 10 possible pairings of the five fitting 
parameters. In each case a range of model solutions {xl < 2.0) for each pair of parameters 
is plotted, while the remaining three parameters are held fixed at their optimum values 
(Table 7). 



